Genome Biology and Evolution

Animal mitochondrial rRNAs are commonly viewed as structurally reduced, yet sponge mt SSU rRNAs range from compact to highly expanded structures. Using nine conserved structural anchors, we compared 216 taxonomically resolved records from four classes and 22 orders, including 16 freshwater Spongillida and 200 marine sponges. Twelve homologous hypervariable substructures were coded as structural types, and their ordered combinations as composite types. We identified 38 structural types and 62 composite types across molecules ranging from 853 to 2,019 nt. Hexactinellida and freshwater Spongillida were each uniform for a distinct composite type but differed markedly in overall structure: hexactinellid mt SSU rRNAs were compact, whereas those of Spongillida were long and contained four to five candidate insertion regions. These results show that a conserved scaffold can accommodate extensive lineage-associated structural variation and provide a practical framework for comparing highly divergent mitochondrial rRNAs.

Chromosome evolution varies widely across Aphidomorpha: several aphid lineages exhibit extensive interautosomal reshuffling, whereas comparisons involving grape phylloxera and Adelges suggest greater chromosome conservation outside Aphididae. However, all published adelgid genomes represent the genus Adelges, leaving conservation across deeper adelgid divergences unresolved. Here, we report a chromosome-scale genome for the pine bark adelgid, Pineus strobi, providing the first genome for the genus Pineus and extending genomic sampling to an early-diverging adelgid lineage. Synteny analyses revealed broad conservation of major linkage groups between Pineus and Adelges across approximately 90 million years. The five grape phylloxera chromosomes also showed broad correspondence to the ten adelgid chromosomes, consistent with a limited number of chromosomal fusions or fissions and relatively little exchange among major linkage groups. These findings strengthen evidence that such extensive interautosomal reshuffling is not characteristic of Aphidomorpha as a whole. We also recovered a complete 25.3 kb mitochondrial genome with expanded, repeat-rich noncoding regions and complete circular genomes for the obligate nutritional symbionts "Candidatus Annandia pinicola" and "Candidatus Hartigia pinicola." Comparisons with Pineus similis symbionts revealed conserved coding capacity and genome-wide synteny, supporting conservation of this dual nutritional symbiosis in pine-associated adelgids. Finally, we recovered the first complete Wolbachia genome reported from an adelgid. Because Wolbachia can induce parthenogenesis, its presence in P. strobi raises the possibility that it contributes to the maintenance or reinforcement of parthenogenesis in a species lacking a viable sexual generation. Together, these genomes provide an integrated resource for aphidomorph evolution and symbiosis. Significance statementSeveral aphid lineages exhibit extensive reshuffling among autosomes, whereas previous comparisons with grape phylloxera and Adelges cooleyi suggested greater chromosome conservation outside Aphididae; however, all published adelgid genomes have represented Adelges, leaving the long-term stability of chromosome organization across Adelgidae unresolved. The first chromosome-scale genome from Pineus reveals broad conservation of major linkage groups across approximately 90 million years of adelgid evolution and clear correspondence with the five chromosomes of the grape phylloxera, supporting chromosome conservation outside Aphididae. The accompanying mitochondrial and symbiont genomes provide an integrated resource for studying aphidomorph genome evolution and symbiosis.

Whole-genome duplication (WGD; or polyploidy) has played a major role in the evolution of many lineages however, our understanding of the processes that shape genome evolution following WGD remains incomplete. While polyploidy duplicates the entire genome sequence, it is rediploidisation that establishes independent duplicated genes. Rediploidisation proceeds through suppression of meiotic recombination across polysomic loci thus restoring disomic inheritance, a process that is not synchronised across the genome. Despite its importance, the mechanisms underlying this process remain poorly understood. The slowly evolving genomes of paleopolyploid Acipenseriformes paddlefish and sturgeon provide an invaluable system for investigating this, as rediploidisation was highly asynchronous in these lineages. In both genomes ohnologs tend to segregate into blocks on the chromosomes according to rediploidisation timing, a pattern that suggests links between chromosomal structure and rediploidisation. Here, we analyse a newly-produced duplicate-resolved paddlefish genome assembly and show a strong concordance between genome rearrangement and rediploidisation timing. We also find that topologically associated domain (TAD) boundaries are associated with rediploidisation block boundaries. Together these results indicate that rediploidisation in acipenseriformes occurred through a process of genome rearrangements that was subject to functional constraints imposed by 3D genome architecture. We investigate the evolution of Hox clusters in these lineages, revealing a previously overlooked duplicate HoxC region in paddlefish, and both ancestral and lineage-specific Hox cluster rediploidisation with substantially different timings. These findings highlight a complex evolutionary history following WGD in Acipenseriformes with implications for understanding short-term adaptations to polyploidy as well as longer-term diversification of lineages.

Background: Chromosome-level genome assemblies are increasingly enabling tests of chromosome evolution, conserved synteny, and sex chromosome conservation across diverse animal lineages, including spiders. Results: Here, we present a chromosome-level genome assembly for the African hermit spider, Nephilingis cruentata, a species with extreme female-biased sexual size dimorphism and a cytogenetically inferred XX2 sex chromosome system. The final Hi-C-assisted assembly spans 1.72 Gbp, with 99.5% of bases assigned to 13 pseudochromosomes, a scaffold N50 of 131.6 Mbp, and a BUSCO completeness score of 98.8%. We annotated 20,021 protein-coding genes, and repetitive elements accounted for 42.7% of the genome. Sex-specific whole-genome resequencing identified Chr02 and Chr07 as candidate X chromosomes based on reduced male coverage, consistent with the expected XX2 system. Using comparative whole-genome alignments across existing chromosome-scale spider assemblies, we also show that sex-linked chromosomes retain broad homologous identity across sampled spider lineages but exhibit lower synteny conservation and greater chromosome-length divergence than autosomes. Conclusions: These results suggest that spider sex chromosomes are conserved in homologous identity but more labile in structure, providing a comparative framework for studying sex chromosome conservation and divergence across Araneae.

The eriophyoid mites (Acari: Eriophyoidea) represent an extreme case of genome streamlining, with genomes averaging just 32Mb, among the smallest of all animals. This clade of mites is highly diverse with more than 4000 species. Their morphologies are specialized for feeding on their host plants, including a simplified worm-like body plan with just two pairs of legs and modified mouth parts that can induce the formation of galls or plant deformities during feeding. Their generally strong host affinities make eriophyoid mites appealing for use as biological control agents, although their short generation times and small genomes could facilitate rapid evolution and impact their efficacy in management programs. Here, we sequenced the genome of the biological control mite Floracarus perrepae, producing a highly contiguous genome totalling just 23.1 Mb, among the smallest of all animals. We also assembled another non-eriophyid mite genome from accidental DNA bycatch (69.4Mb). We placed this new genomic resource in a phylogenetic context to reveal that mites have highly dynamic genome evolution, with a significant trend in genome downsizing in the eriophyoids. Our results suggest that this streamlining is associated with non-genic elements such as the suppression or excision of retrotransposons and purging of introns. As new sequencing techniques become available, novel genomic resources for tiny organisms such as F. perrepae will be more readily accessible, facilitating both fundamental genome evolutionary biology and applied sciences such as biological control programs which use eriophyoid mites for the management of invasive species.

The European leaf-toed gecko (Euleptes europaea) is a small, nocturnal gecko endemic to the western Mediterranean. As a phylogenetically distinctive member of the Gondwanan family Sphaerodactylidae, it represents an important species for studying Mediterranean island biogeography, adaptation, and reptile genome evolution. The species also occupies a key position for investigating the evolution of sex chromosomes, as geckos exhibit remarkable diversity and frequent transitions in sex-determination systems. We present a chromosome-level genome assembly of Euleptes europaea generated as part of the Vertebrate Genomes Project. The 1.8 Gb assembly has a scaffold N50 of 102.3 Mb (contig N50 27 Mb), with 21 chromosome-scale scaffolds corresponding to the known karyotype (2n = 42). The primary assembly has a BUSCO completeness of 97.80% (95.60% as single-copy), a k-mer completeness of 96.00%, and a k-mer quality value (QV) of 61.20. Repetitive elements account for 53.20% of the genome and genome annotation identified 18,633 protein-coding genes. This high-quality reference genome will facilitate studies of genome evolution, island adaptation, and sex chromosome evolution across geckos and other reptiles.

We present a high-quality chromosome-level reference genome for the Saimaa ringed seal (Pusa saimensis), an endangered freshwater pinniped endemic to Lake Saimaa, Finland. The assembly spans 2.353 Gb and comprises 15 autosomes together with the X and Y sex chromosomes. Using Oxford Nanopore Technologies (ONT) long-read sequencing and Hi-C scaffolding, we achieved a telomere-to-telomere assembly for all chromosomes, except the Y chromosome. Genome annotation identified approximately 21,800 protein-coding genes, consistent with other mammalian genomes. Assembly completeness was high, with BUSCO analysis recovering 99.6% of expected complete single-copy genes (98.2% single-copy and 1.3% duplicated). Comparative analyses revealed a highly conserved chromosomal architecture, with only minor syntenic differences relative to other pinniped chromosome-level assemblies. Previously described cytogenetic fusion events in Phocidae were confirmed (chromosomes 2 and 7). A translocation between chromosomes 6 and 7 distinguishes phocids from the otariids. In general, more distantly related taxa exhibit an increasing degree of intrachromosomal rearrangements. Notably, we identified a large intrachromosomal rearrangement on chromosome 2 that appears specific to the Saimaa ringed seal. Phylogenomic analysis based on 9,226 single-copy orthologues placed the Saimaa ringed seal as a sister lineage to the Baltic ringed seal (Pusa hispida botnica), while confirming also other established evolutionary relationships among pinnipeds. Comparative gene family analysis between the Saimaa ringed seal and the closely related grey seal (Halichoerus grypus) revealed lineage-specific differences driven by a limited number of gene families. In the Saimaa ringed seal, expansions were observed in ion transport, cytoskeleton, and regulatory genes, potentially reflecting adaptation to freshwater conditions. In contrast, the grey seal showed expansions in olfaction, immune-and spermatogenesis-associated gene families, including MAGE/MIA genes, consistent with differences in ecology and mating systems. This reference genome provides an important resource for studies of pinniped genome evolution, as well as conservation and population genomics of the Saimaa ringed seal, facilitating future work on genetic diversity, inbreeding, mutational load and adaptive potential in this highly endangered species.

Aspergillus species are ecologically diverse and deeply entangled with human health and industry. A. fumigatus and A. flavus are the two principal species of invasive aspergillosis. A. niger and A. oryzae, on the other hand, are responsible for global enzyme production, organic acid production, and koji-based fermentation industries. The question of whether these similar phenotypes share the same genomic mechanisms across the genus is not yet understood. To address this, we constructed per-species pangenomes for the four Aspergillus species (929 initial genomes filtered to 210 ANI-verified, high-quality assemblies for a total of 88 A. fumigatus, 70 A. flavus, 33 A. oryzae, and 19 A. niger assemblies) alongside a genus-level pangenome of 15,163 orthogroups, and conducted phenotype-labeled pan-genome-wide association studies (pan-GWAS) with kinship correction across all species. Pan-GWAS identified up to 117 significant orthogroup presence/absence associations per species-phenotype comparison. However, convergence analysis showed that among the 92 and 62 distinct gene families significant for human pathogenicity in A. fumigatus and A. flavus respectively, the two species seldom agreed on whether the pathogenicity was associated with the enrichment or the depletion of a specific gene family. Convergence analysis of the functional annotations also yielded zero significant results at FDR < 0.05. A literature-curated gene panel analysis also showed that a species labeled pathogenic and another labeled GRAS carried the same aflatoxin and virulence genes, suggesting that gene presence alone cannot readily explain their phenotypic differences. Instead, we propose that niche adaptation operates through the use of the pangenomic rare genome. Reclassifying rare genes by homology identified truly rare subsets (156 to 391 orthogroups per species) distinct from paralogs and gene fragments. Human-pathogenic strains showed significant rare genome expansion of 2.44-fold for both A. fumigatus and A. flavus (kinship corrected, p = 6.6 e-08). Conversely, industrial strains showed rare genome contraction where both A. niger and A. oryzae industrial strains carried 0.57-fold (kinship corrected, p = 0.015) fewer rare genes than their non-industrial counterparts. Hence, we claim that Aspergillus niche evolution proceeds through directional rare genome changes, where there is expansion under pathogenic selection, and contraction under industrial domestication. The rare genome, often discarded as noise, may represent the primary evolutionary source for clinical and biotechnological adaptation in this genus.

Chloroplast or plastid genomes are essential resources for studying the evolution and diversity of algae and land plants. Although thousands of plastid genomes have been sequenced, their full potential has not been realised; derived resources such as orthogroup databases and reference datasets for metagenomic profiling remain underdeveloped. We present the ChlORIS database to address these problems across all algal phyla. From 2,254 publicly available algal plastid genomes, after dereplication we clustered 2,531 orthogroups from the annotated proteins and selected 496 orthogroups with consistent gene naming, enabling cross-genome comparisons of homologous plastid proteins. We further selected 224 core orthogroups, each containing more than 10 protein sequences, for which we produced score-calibrated hidden Markov models (HMMs), multiple sequence alignments and predicted protein structures. The value of these resources for phylogenomics is demonstrated through a large-scale plastid phylogeny of 859 taxa spanning all major algal lineages. We characterised the protein HMMs by cross-referencing them to Pfam domains and calibrated score cutoffs for reliable detection. The metagenomic database, HMM library, nucleotide and amino acid alignments, predicted structures and protein metadata, cross-linked to UniProt and InterPro (Pfam), are openly available on the ChlORIS website at https://chloris.codeberg.page/.

Transposable elements (TEs) constitute a large fraction of mammalian genomes yet their contribution to variation among individuals within natural populations remains largely unexplored. While most TE insertions are deleterious, some may be beneficial and contribute to adaptation. We characterized TE variation and assessed its potential adaptive role using long-read whole-genome sequencing of wild-caught house mice (Mus musculus domesticus) sampled from two populations inhabiting contrasting temperate and tropical environments and differing in morphology, physiology, and behavior. We sequenced 10 mice from each population and created highly contiguous de-novo genome assemblies for each individual, allowing us to identify TEs that are not present in the mouse reference genome and to characterize individual variation. By performing manual TE curation, we identified 506 non-redundant TE consensus sequences among all mice. On average, each wild mouse genome contained 1.47 million TE insertions, ~4% of which were polymorphic among individuals. A small fraction of these polymorphic TE insertions were present in high frequency in just one of the populations, consistent with positive natural selection. Using liver RNA-seq in natural populations and in laboratory crosses, we studied gene expression at genes adjacent to polymorphic TEs. This identified a small set of TEs that are associated with the expression of nearby genes in a population-specific manner, nearly all of which showed independent signatures of positive selection. Together, these results provide the first detailed assessment of TE variation in natural populations of house mice and identify a small set of TE insertions that likely contribute to environmental adaptation.

Alternative splicing (AS) can generate multiple RNA isoforms from a single gene and is thought to contribute to phenotypic divergence, including differences between the sexes. Studies in several organisms have documented sex differences in splicing, however these have been largely reliant on short-read RNA sequencing which requires complex algorithms to assemble full-length transcripts and may underestimate both isoform diversity and sex differences in splicing. We used long-read, single molecule RNA-Seq to build a more complete catalog of sex-biased splicing in Poecilia reticulata, a focal species for studies of sexual dimorphism. Pairing long-read sequencing with isoform-level analyses, we identified a sixfold higher proportion of sex-biased splicing genes (37%) compared with short and long-read event-level approaches (6%, 11%). AS was common (70% genes) but only 54% of isoforms produced unique open-reading frames (ORFs). We found that males exhibited greater isoform richness than females in both tail and gonad tissues but produced a smaller proportion of isoforms with unique ORFs, suggesting that much of the increased isoform variation is unlikely to expand proteomic complexity and may instead reflect stochasticity during splicing rather than intentional transcriptional products intended for translation. Despite widespread AS, we found only 2.3% of genes exhibited sex-biased isoform switching, and only 52% of these switches generated distinct sex-biased ORFs. Together, our long-read data suggest that although isoform diversity is more extensive than previously appreciated, most alternative isoforms are unlikely to generate novel proteins. Instead, a relatively small number of sex-biased isoforms may disproportionately contribute to proteomic divergence between the sexes.

The Mutational Hazard Hypothesis (MHH) predicts that reduced effective population size weakens purifying selection and promotes transposable element (TE) accumulation. Because effective population size is difficult to estimate across broad taxonomic scales, body mass, generation time, and dN/dS are often used as proxies. We analyzed TE landscapes across 167 sauropsid genomes to test whether these proxies predict genomic TE proportion consistently across phylogenetic scales. All three showed strong scale dependence. Body mass was positively associated with TE proportion across clades, but this relationship broke down within clades: only snakes retained a significant negative Pearson correlation after multiple-testing correction, whereas the corresponding phylogenetically corrected slopes were not significant. Generation time showed a strong pooled association that disappeared within every major clade in both Pearson and PGLS analyses. The pooled dN/dS-TE relationship also disappeared after accounting for body mass and generation time in a structural equation model, with a robust within-clade association retained only in turtles. Lineage-specific TE dynamics, including LTR expansion in sea snakes, were not captured by these proxies. These results show that commonly used MHH proxies mainly reflect clade-level structure rather than consistent within-lineage mechanisms.

Butternut (Juglans cinerea) is a North American hardwood in steep decline, driven largely by an introduced fungal pathogen that causes butternut canker disease. Tolerance exists within the genus, with the Asian Juglans ailantifolia comparatively tolerant and hybridizing readily with J. cinerea, yet the genomic basis of that difference is poorly resolved. Chromosome-scale reference genomes are presented for J. ailantifolia and for a second J. cinerea accession from the species primary range, with independent scaffolding and uniform annotation across all compared assemblies. Gene family evolution was assessed across twelve Juglandaceae genomes, one per species, and a four accession pangenome was built from J. ailantifolia and J. mandshurica representing section Cardiocaryon and two J. cinerea accessions, Brenton and New Brunswick, representing section Trachycaryon. Neither J. cinerea accession is a putative tolerant selection, so their contrast provides a null against which section-level differences are calibrated. Spanning two sections separated by approximately 30 My, the pangenome partitions many tandem arrays by accession instead of by gene, so copy number is interpreted after merging homology nodes representing the same gene. Section-level variation resides primarily in gene presence and absence, where pathogenesis-related proteins, cysteine-rich receptor-like kinases, wall-associated kinases and lectin-domain receptor kinases are over-represented in the lineage-specific complement of both sections, indicating rapid bidirectional turnover. A smaller set of families is asymmetric, with chitinases and chitin-binding proteins and NB-LRR receptors biased toward the tolerant section and dehydrins, pectin-modifying enzymes and the CBF regulon toward the susceptible one. The largest copy number difference in the dataset separates the two conspecific accessions, at a senescence-associated cysteine protease. The two butternuts are further distinguished by their demographic histories and by the functional composition of their accession-specific gene content: the northern New Brunswick accession retains lower heterozygosity, a distinct coalescent trajectory, and enrichment for calcium transport across channels, pumps and homeostasis, for salt stress response regulation, and for raffinose family oligosaccharide biosynthesis, spanning both the perception and the cryoprotective effectors of freezing tolerance. These findings support recognition of the northern population as a cryptic glacial refugium and a conservation priority.

High-quality reference genomes provide an essential foundation for elucidating the molecular basis of organismal ecophysiology. Here, we sequenced and assembled chromosome-scale genomes of two Heterosigma akashiwo strains isolated from coastal waters of Japan and France. The assembly sizes were 1.18 Gb and 1.43 Gb for the Japanese and French strains, respectively. The scaffold N50 of the Japanese strain assembly was 66 Mb, whereas the one of the unscaffolded French strain assembly was 33 Mb. To our knowledge, these assemblies represent among the largest and most contiguous genome resources currently available for members of the Stramenopiles (Ochrophyta). Evidence-based gene prediction in the Japanese strain recovered approximately 90% of conserved stramenopile core genes, indicating a highly complete gene repertoire, and was complemented by extensive functional annotation. In the French strain, homology-based gene prediction recovered approximately 80% of conserved core genes. Comparative genome analysis revealed extensive synteny conservation between the two strains, although several putative duplication and translocation events were detected. These genomic resources provide a robust framework for investigating the molecular, cellular, and ecological mechanisms underlying the physiology, adaptation, and bloom-forming capacity of H. akashiwo.

Archaic introgression has shaped the evolutionary history of Eurasian populations, yet Central Eurasian region remains understudied despite being at the crossroads of ancient human migration. Here, we analyzed the whole-genome data of five Central Eurasian (CE) individuals from Early Bronze Age (EBA) and five present-day CE individuals to characterize the archaic introgression landscape. We estimated that archaic introgression from Neanderthal and Denisovan archaic hominins comprises approximately 2.2% of the Central Eurasian genomes. Both amount and chromosomal distribution of archaic introgression remained largely unchanged between the EBA and present-day CE individuals. Putative introgressed fragments matching the Altai Neanderthal and the Altai Denisovan were retrieved. Our results suggest that while the archaic introgression levels seemingly remained stable over the past several thousand years, larger modern CE genomes panels will be required to fully characterize the genomic landscape of archaic ancestry in the region.

Protein evolution is a walk in the evolutionary space directed by mutation and selection. While functional and structural constraints serve as the main determinant of amino acid substitution in most proteins, synthesis cost and mutational bias can also alter the direction and rate of amino acid evolution, especially in systems experiencing relaxed selection. Here, we focused on the highly expressed plastid ribosomal proteins (PRP), which comprise 58 conserved proteins encoded by both plastid and nuclear genomes. Relaxed selection has been repeatedly identified in three distantly related plant lineages, providing a valuable comparative framework to investigate the significance of synthesis cost and mutation. We first demonstrated that the hemiparasitic tribe Cymbarieae (Orobanchaceae) represented a new case where concerted cyto-nuclear rate elevation occurs in their PRP. Further investigation revealed convergent shifts in amino acid composition in all four plant lineages attributable to arginine-to-lysine and methionine-to-isoleucine/valine/leucine substitutions. The replacement residues were biophysically similar but had lower molecular weight and shorter side chains, which significantly destabilized protein folding as demonstrated by protein structure modeling. We found that the composition shifts ran counter to the expectation of mutational bias but were consistent with the expectation of synthesis cost minimization, which is potentially adaptive for highly expressed PRP. Further, cost minimization significantly influenced all conservative substitutions between biophysically similar amino acids but was absent in non-conservative substitutions. We thus propose cost minimization as a secondary selective drive for protein evolution in PRP, unmasked in lineages and sites with relaxed selection on their function.

Genetics may help address the biodiversity crisis by providing information about genetic diversity and temporal changes in demography for species of interest. Advances in whole-genome sequencing create new opportunities for demographic analysis, even based on the two copies of a genome found in a single diploid individual. The Vertebrate Genomes Project (VGP) is generating high-quality, chromosome-level reference genomes across the full range of extant vertebrate species, with its first phase delivering assemblies spanning approximately 95% of vertebrate orders. Using 512 diploid VGP genomes, we quantified intra-species heterozygosity, runs of homozygosity (ROH), and inferred past effective population sizes (Ne) with the Pairwise Sequentially Markovian Coalescent (PSMC). Threatened species are more likely to exhibit lower heterozygosity and longer ROH, though there is large variation in both measures across all IUCN categories. Interestingly, PSMC suggests that estimated historical Ne several thousand generations ago is a better predictor of threatened status than the present day estimate. Co-analysing with life history traits, we found that marine species tend to have lower ROH content, while fossorial species show significantly higher inbreeding levels. Indeed, habitat and foraging strata are much stronger predictors of IUCN status than genetics, with estimated historical Ne providing a small but significant amount of additional information. Together, these results suggest that, while measures of genetic diversity are correlated with IUCN status, much of that correlation may derive from ecological factors such as habitat, with only a relatively small direct contribution. Nevertheless, reference genomes like those generated by the VGP can yield valuable information, like historical Ne, while facilitating population monitoring and management for species of interest.

Insect-parasitic nematodes are widespread and often significantly reduce host fitness, yet we know surprisingly little about most species. Contortylenchus reversus is a hemocoel-inhabiting parasitic nematode that infects Dendroctonus bark beetles, notably impacting host mobility and fecundity. We first detail a chromosome-scale genome assembly of C. reversus, recovered serendipitously from a sequencing project targeting a host (Dendroctonus ponderosae). We assembled the 79.4 Mb genome into nine linkage groups and, through transcriptome-aided annotation, identified 11,244 protein-coding genes. Synteny comparisons with the only relatives for which there are complete assemblies reveal extensive chromosomal rearrangements and extreme loss of gene collinearity suggesting these insect-parasitic nematodes may have exceptionally malleable genomes. Using this genome assembly and repurposed reduced-representation genomic data from 707 D. ponderosae individuals, we investigated infection frequencies and population structure, identifying infection rates ranging from 0% to 50% across 18 geographically widespread collection sites. Population structure of C. reversus appears broadly concordant with the structure of the host beetle, suggesting a shared evolutionary history, while genetic variation (nucleotide diversity) in the parasitic nematode is highly reduced in comparison to its host. These results offer insights into its population genetics, host associations, and the evolutionary dynamics of a nematode-beetle interaction and highlight how genomic bycatch can reveal previously hidden details about an important species in a complex community.

Chromosomal inversions are increasingly recognized as important drivers of local adaptation and ecological divergence because they suppress recombination and maintain adaptive allele combinations despite ongoing gene flow. However, the eco-evolutionary conditions favouring the establishment of such indirectly adaptive inversions, as well as the genomic features that distinguish them from other inversions remain incompletely understood. In this study, we investigated these questions in a wild sunflower system comprising two species: Helianthus debilis and Helianthus praecox, which exhibit diverse ecotypes and varying degrees of geographic overlap across Texas and Florida in the USA. To resolve the evolutionary relationships between and within these species, and identify potentially adaptive inversions, we generated haplotype-resolved reference assemblies and integrated comparative and population genomic analyses. We identified three major genetic clusters that only partially corresponded to taxonomic classifications. We further detected 156 inversions across the genome, 11 of which showed signatures consistent with a role in local adaptation. Notably, nine of the 11 putatively adaptive inversions were found in sympatric Texas populations. Together with a similar enrichment of inversions in genome assemblies from sympatric versus allopatric populations, our results suggest that inversions are more likely to evolve in heterogeneous environments with ongoing gene flow than in allopatry. Lastly, locally adaptive inversions were generally larger, contained more genes, and showed greater sequence divergence between haplotypes than other types of inversions. Our findings provide empirical support for the role of gene flow in promoting inversion establishment and identify genomic characteristics associated with indirectly adaptive inversions.

Transposable elements (TEs) are ubiquitous repetitive DNA sequences that can mobilise within genomes and may modulate gene expression in an environment-dependent manner. TEs and the safeguarding epigenetic machinery targeting them, can be tuned by environmental fluctuations to influence gene expression by inducing genomic, epigenetic, and transcriptomic changes. Yet, the degree to which TE-driven molecular diversity translate into inter-individual phenotypic variation vs accumulating without any phenotypic consequences remains unclear. Here, we used five populations of genetically engineered Drosophila melanogaster flies that carry variable TE content but share an otherwise identical genetic background to test the phenotypic consequences of the early stages of TE accumulation. Phenotypic screenings across 17 traits (fertility-related traits, life-history traits and stress resistance tests) revealed significant differences between the populations (e.g. reduced hatchability). We also observed a notable increase in intra-population phenotypic variation for the heavily TE-burdened populations across a wide panel of traits. These results suggest considerable TE-driven inter- and intra-population phenotypic variation. Further investigation revealed that variable TE contents can influence the response to environmental changes, positioning TEs as drivers of environmentally-induced phenotypic variation in a system deprived of other sources of genetic variation. These results provide empirical evidence that TEs contribute to the heterogeneity of the environmental response and therefore represent an underlying mechanism of phenotypic variation.